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ABSTRACT 

Helium core white dwarfs (WDs) with mass M < 0.20 Mq undergo several Gyrs of stable hydrogen 
burning as they evolve. We show that in a certain range of WD and hydrogen envelope masses, these 
WDs may exhibit g-mode pulsations similar to their passively cooling, more massive carbon/oxygen 
core counterparts, the ZZ Cetis. Our models with stably burning hydrogen envelopes on helium cores 
yield g-mode periods and period spacings longer than the canonical ZZ Cetis by nearly a factor of 
2. We show that core composition and structure can be probed using seismology since the g-mode 
eigenfunctions predominantly reside in the helium core. Though we have not carried out a fully 
nonadiabatic stability analysis, the scaling of the thermal time in the convective zone with surface 
gravity highlights several low-mass helium WDs that should be observed in search of pulsations: 
NLTT 11748, SDSS J0822-f2753, and the companion to PSR J1012+5307. Seismological studies of 
these He core WDs may prove especially fruitful, as their luminosity is related (via stable hydrogen 
burning) to the hydrogen envelope mass, which eliminates one model parameter. 

Subject headings: stars: white dwarfs — stars: oscillations 



1. INTRODUCTION 

White dwarfs (WDs) are observed to pulsate in nor- 
mal modes of oscillation (g-modes) which are determined 
by the structure of the stellar interior and atmosphere 
(jWinget k Keplej [20081) . Those with hydrogen atmo- 
spheres exhibit pulsations when they enter the ZZ Ceti 
variable (DAV) instability region, a discrete strip in the 
Teff-logg plane that spans 11,000K < T^s < 12,250K 
at \ogg w 8.0. The ZZ C eti strip has been inves- 
tigated both th eoretic ally (Bra ssard fc Fontaine! 119971 : 
iWu fc GoldreichI 119991: iFbn tainc ct al. 2003) an d em- 
iricallv (iWesemael etal.1 11991,: .Mukadam efaLl 120041: 



pir: 



iCastanheira et al.ll2007l : iGianninas et al.ll2007| ). To date, 
all known ZZ Ceti pulsators have masses > 0.5 M^g, im- 
plying cores composed of carbon, oxygen, and heavier 
elements. 

Lower mass (Af < 0.5 Mq) WDs with nearly pure he- 
lium cores are made on the red giant branch (RGB) when 
core growth is truncated before reaching «0.45— O.47M0 
(logg « 7. 67 at T^s » 11, 500 K), prior to the helium 
core flash (ID'C ruz et all 119961: IDominguez et all 119991: 
iPietrinferni et al. 2004 ). Two modes of envelope mass 
loss can cause this: strong winds or binary interaction. 
Significant mass loss due to stellar winds in high metal- 
licity systems may strip the H envelope, preveii ting the 
He core flash (ID'Cruz et aLlfloM IHansenl [2005h . Com- 



mon envelopes induced by binary interactions also lead 
to sig nificant mass loss (jlben fc LivicillT993l : iMarsh et al.l 
11 9951 ) ■ and make very low-mass He WDs (M < 0.2 M©) 
when the binary i nteraction occ u rs at the base of the 
RGB ( van Kerkwijk.eLalJ 119961 : iCallanan etHI 119981 : 
iBassa et al.l 120061) . Thus, He is the expected core com- 
position for WDs below «0.45— O.47M0. However, 
very little direct evidence exists of He cores. The 
overbrightness of ol d WDs (iHanseiJ 120051) in the star 
cluster NGC 6791 (jBedin et all 120051 ) presents possi- 
ble evidence. The detection of low \ogg young WDs 
make a plausible a rgument for the old WDs to be He 
core (jKalirai et al.l I2007D. however, other possible ex- 
planat io ns remain (iDelove fc Bildstenll2002l : iBedin et al.l 
l2008allbl : iGarcia-Berro et al.ll2010l ). 

Asteroseismology offers the possibility of directly con- 
straining the He core composition in these low-mass 
WDs, as the g-mode periods provide information on 
WD mass, mass of H envelope, and core composi- 



tion jCorsico fc Benvenutol 120021 : ICastanheira fc Kepleii 
120081) . Theoretically, these principles have been a pplied 
to C/0 versus 0/Ne core WDs b ylCorsicoet all (1200 " 
and O versus He core WDs by lAlthaus et a l.' f2004' 
However, to carry this out, we need to find pulsating 
He core WDs, something that has yet to occur. 

Past studies have illuminated a dichotomy i n the 
evolution of the He core WDs (jDriebe et al.l 119991 : 
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iSerenelli et al.ll2002HPanei et al.l2007D that impacts their 
seismic properties and prevalence as pulsators. For 
masses >O.2M0 (dependent upon metaUicity) the H en- 
velope experiences a multitude of H shell flashes that 
reduces its mass, eventually allowing the WD to cool 
rapidly. Such objects traverse the extrapolated ZZ Ceti 
instability strip in ~ 10-100 Myr (jPanei et al.ll2007[) . al- 
lowing for an investigation of their H layer mass, and 
confirmation of pure helium core. However, there are 
presently no known WDs in, or near, the extrapolated 
strip for masses in th e 0.2 M© < M < 0.5 Mq range 
(jSteinfadt et al.l 120081 ). inhibiting such research. Less 
massive (< 0.2 M©) helium WDs have a different evo- 
lution, undergoing stable H burning for Gyr, slowing 
their evolution to rates that may yield more in the ex- 
trapolated instability strip. However, the presence of a 
thick, actively burning hydrogen layer requires new seis- 
mic modeling, especially for the eventual assessment of 
the of the instability strip for these u nusual WDs. The 
recent discovery of thr ee such objects (jKawka fc Venned 
l2009t iKilic et alll2010[ ) makes our work quite timely. 

Motivated by a desire to study the pulsational prop- 
erties of these long-lived systems, we begin in Sec- 
tior(2] by constructing He core WD models with sta- 
ble H burning shells fiexible enough for seismic inves- 
tigations, and compare to the results froin evol ution- 
ary codes (|Serenelli et al.ll2002t iPanei et al.ll2007[ ). We 
discuss the unusual seismic properties of these objects 
in Section [Sj where we calculate their adiabatic mode 
structure and, using an approxim ation for the in stabil- 
ity criterion of Brickhill's theory ()Brickhil]||1991[ ). high- 
light the potential location of the He core WD insta- 
bility strip. In Section |4l we suggest a few intrigu- 
ing pulsation candidates amongst the very lowest mass 
WDs (van Kerkwiik et al. 1996; Ca ll anan et al. 1998; 
Bassa ct all 120061: iKawka et al.l l2006t iKawka fc VennesI 
2009; .Kilic et al.ll2010D . where our simple models apply. 



We close by highlighting the need for future work, espe- 
cially if observations of our suggested targets yield the 
first pulsating, low-mass. He WD. 

2. HYDROGEN BURNING MODELS 

The < 0.2 Mq He core WDs of interest for our work 
are undergoing stable H burning via the PP chains in a 
low-mass shell. The solid/dash e d line s in Figure [1] are 
the models from ISerenelli et al.l (|2002[ ) and IPanei et al.l 
([2007). These lines transition to solid when 90% of the 
luminosity is generated from the PP chains. The solid 
circles are the locations of the observed WDs of interest 
here, and are clearly in a region dominated by PP-chain 
burning. For this reason, we construct models with only 
PP-chain burning. In addition, at these late times, the 
diffusive timescale at the burning zone is much shorter 
than the age and the WD core temperature is set by 
that in the stable burning layer. For these reasons, the 
prior evolution of a stably burning WD does not affect 
its properties at this stage of evolution. 

We construct models of a stably burning H envelope 
on an He core, by solving the equations of hydrostatic 
balance, heat transport, energy generation, and mass 
conservation. Between the H and He layers the most 
important physics in our models is the chemical profile. 
In these layers diffusive equilibrium is valid as the evolu- 
tionary timescale (dominated by nuclear burning, ~ 1- 
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against those of IPanei et al.l IJ2007I ) . 0.16 M©, and ISerenelli et al.l 
(iSni), 0.20, 0.25, 0.35, and 0.45 Mq, the thin black lines. 
Our tracks de viate significantly at low gravity from those of 
ISerenelli et al.l ((2002). The solid circles are several candi- 
date < 0.20 Mpi WD pulsation candidates llvan Kerkwiik ct ^ 
TgoS; iCallanan ct al. 1998; Bassa et aLl [200g: IKawka et al.l 1200^ 
Kilic ct al. 2010; Kawka & Vcnncs 2003). Note their location on 
the higher gravity, past the T(,ff turn-around portion of the models. 
The dashed lines in the ScrcncUi ct al. (2002) models denote where 
the CNO luminosity is greater than 10% of the PP chains. This 
explains the deviation as our models only have PP-chain luminos- 
ity, appropriate to the most relevant regions at higher gravity. Five 
Gyr have elapsed between the triangle and square marks and the 
square and circle marks. The gray area approximates the location 
of the C/O WDs with the vertical dashed lines denoting the empir- 
ical ZZ Ceti instability strip (Gianninas ct al. 20071). The dotted 
line above the empirical instability strip denotes the gravit y of the 
lowes t known gravity ZZ Ceti, HE 0031-5525 I^Castanheira et al.l 
[200^) . 



10 Gyr) is significantly longer than the diffusive timescale 
(^ 10- 100 Myr) over a pressure scale height. We derive 
the equilibrium electric field by assuming each species is 
in hydrostatic balance wit h gravity and the electric f orce 
and charge neutrality (see iChang fc Bildstenl (|2003f l for 
more detailed derivation). With this electric field we 
generate an additional differential equation for one of 
the chemical species (charge neutrality gives the rest) 
to be simultaneously solved with the equations of stel- 
lar structure. Given our set of differential equations and 
boundary conditions, our model reduces to two parame- 
ters, total mass and total H mass, fewer than those for 
passively cooling C/O WDs which require total mass, H 
mass. He mass and surface temperature. 

Our models contain zero metallicity. Addi- 

tional elements require additional differential equa- 
tions for their diffusive profile for which equi- 
librium conditions may not exist. Therefore, 
the PP-chains {p{p,e+Uef}i{p,-f)^}le{^}le,pp)'^}le 
or ^lle{'^Re,-/yBe{e-,Pe{l)yU{p,ayRe or 
^Be(p,7)*B(e+i^e)^Be 2 "^He) are the only source 
of nuclear energy. We assume ^He has reached its 
equilibrium abundance, peaking at ~ 10""' to 10~^ by 
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mass. We generate evolutionary tracks by conserving 
total mass; nuclear burning converts envelope mass 
into core mass. Figure [1] com pares our mode l s with 
the time-dependent m odels of iSerenelli et al.l ()2002l ) 
and Pane l et al.l (|2007[ ). exhibiting discrepancy at low 
gravity (large envelope mass) but excellent agreement 
at high gravity (low envelope mass). This is attributed 
to large CNO luminosities at high envelope mass in the 
non-zero metallicity models of ISerenelli et al.l ()2002fl . At 
low envelope masses, the core temperatures have lowered 
to w 10^ K and CNO elements h ave diffused out of the 
burning region (jPanei et al.ll2007l ). Therefore, PP-chain 
luminosity dominates so we expect our models to be 
valid in this regime, where, as we show, the candidate 
He core pulsating WDs are likely to be found. As 
nuclear burning determines the evolution of our models, 
the timescales are long, of order several Gyr. Figure [T] 
illustrates the evolution of our models over 5 (triangles 
to squares) and 10 Gyr (triangles to circles) from the 
point of maximum Toff. 

The microphysics, opacities, equation of state, and nu- 
clear energy generation, are all handled by the Modules 
for Experiments in Stellar Astrophysics (MESAQ code, 
developed by B. Paxton et al. (2010, in preparation). 
Within MESA, the opac it ies ar e dra wn from OPAL 



J 11 



( Inglesias fc RogersI 119931 I1996D . the iFerguson et"al 
( 2005") low-temperature tables, and the iCassisi et al 
(2007) electron conduction table s. The equation 
of state is derived from OPAL (Rogers & Navfono" 
I2002f ). low-temperature SCVH (Saumon et al. 19951 
and fully ionized high te mperature and density HELM 
([Timmes fc Swestvl 120001 ). Nuclear energy gene ration is 
calcul ated using the techniques developed by iTimmesI 
(p99h . 

3. NON-RADIAL PULSATION ANALYSIS 

To analyze the non-radial pulsational properties of our 
He-core WD models we perturb and linearize the fluid 
equations of momentum, energy, and mass conservation. 
We set the transverse wavenumber of order (. as kf^ = 

3.1. WKB Approximation 

The star is divided into regions of wave propaga- 
tion and evanescence. In the propagation zone, the 
wavelength k~^ is much smaller than the characteris- 
tic length scales associated with the background, such 
as the radius r, near the center, and the pressure scale 
height \p = p/{pg) near the surface. This allows for 
the WKB approximation where all state variables are 

ocexp(j/ drkr). Neglecting perturbations on the grav- 
itational field (the Cowling approximation) we further 
reduce the linearized pulsation equations into the dis- 
persion relation. 



r^r — 



N^ 



(1) 



where uj is the frequency of pulsation, N is the Brunt- 
Vaisala frequency, and Cs is the adiabatic sound speed. 
For propagating waves. Equation ([T]) defines the reso- 
nant cavity for waves of two types. When ui > N and 

^ http://incsa.sourceforge.net 
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Fig. 2. — Propagation diagram (upper panel) for our MwD 



0.17 M0 and Me; 



3.15 X 10 •^ Mq model. Gravity mode (g- 



mode) pulsations exist below both the Brunt- Vaisala frequency 
(solid curve) and the Lamb frequency (£ = 1; dashed curve). The 
thin horizontal lines delineate the locations of the n = 1, 5, and 
10 modes for £ = 1. Note the large bump in the Brunt-Vaisala 
frequency due to the Ledoux contribution (Equation |4]l that aligns 
with the composition transition region (bottom panel). The middle 
three panels show the eigenfunction solutions of §;, (dashed curves) 
and mode energy (solid curves) for the n = 1, 5, and 10 modes 
for 1 = 1. The gray vertical lines denote the boundaries of the 
WKB propagation cavity. These show the bulk of each mode to 
reside in the core, below the H/He transition, with only a couple 
nodes existing in the envel ope. This mo del is our closest fit to the 
properties of NLTT 11748 IIKawka fc Venncs 2009) . 



Cskfi (the Lamb frequency), waves propagate as sound 
waves (p-modes) where pressure provides the restoring 
force. When uj < N and Cskf^ waves propagate as grav- 
ity waves (g-modcs) where gravity provides the restoring 
force. [Brickhill (1 985 ) showed that convective driving 
could drive the amplitude of 5-mode pulsations to an ob- 
servable l evel; these are the pulsati ons observed in the 
ZZ Cetis (jWarner fc Robinson|[T97^ . 

We approximate the frequencies of propagating 
g-modes using the WKB quantization condition 
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J °"' drkr = rnr. Under the assumptions uj <^ N and 
w ^ Cskh, Equation ([T]) gives, 



UJn, I 



^W+T) 



nn 



Taut ^^ 

—N, 
r 



(2) 



where the integral is bounded by the frequency depen- 
dent resonant cavity. Here rin and Tout are the radii 
where uj = N and uj = Cgkh, respectively; see Figure [5] 
for illustration. Under these assumptions, the derived 
mode periods are only accurate for large radial order, 
n ^ 1 . Figure [2] shows a propagation diagram for an 
MwD = 0.17Mq and Me„v = 3.15 x 10~^ Mq model. 
As is evident, the resonant cavity for the higher order g- 
modes samples much of the core and envelope while the 
lower orders are most affected by the transition region. 
It is obvious that the contribution of the transition re- 
gion is quite important, therefore, close attention must 
be paid to the Brunt-Vaisala frequency. For changing 
composition, the Brunt-Vaisala frequency is. 



where, Xp = 

/-I 



— (Vad - V) 
^P 



B 



(3) 



and, 
d\x\p 



OlnP/ainplj, _j^^j,Xt 
dX, 



dX, 



T^P,{X,^,} 



dlnP' 



(4) 



is the compositionally dep endent Ledoux term (modified 
from iBrassard et al.lTl991[ ) which accounts for the bulk 
of the bump in the Brunt-Vaisala frequency at the com- 
position transition zone in Figure [5J 

3.2. Numerical Analysis 

To obtain more accurate mode periods that are not 
restricted to the high radial orders as our WKB analy- 
sis is, we must turn to the boundary value problem for 
adiabatic pulsations. The eq uations and method are de- 
scribed in lUnno et al.l |1989). Perturbation of the equa- 
tions reduce to the following three equations for the pres- 
sure perturbation Sp = pijj, the radial Lagrangian dis- 
placement ^r (the transverse Lagrangian displacement is 
£,h — i'r/uP'), and the gravitational potential perturba- 
tion 56: 



dip , 2 



d£.r 
dr 



N' 



dS(t> 



N')(r + —i^~-J^ 



1 d 



, dScjj 
dr 



1 



V' + -^H, 



(5) 
(6) 



N' 



= AttGp\ —tp + C 



klM. (7) 



At the center of the star, requiring the variables to be 
finite leads to the scalings ip oc r^ and ^,. ex r^~^. Placing 
these scalings into Equations ([5]), ([6|), and ([7]) leads to 
the central boundary conditions. 



UJ^^r^e 



ll^ + i 



dr 



(8) 
(9) 



-(^+1)^ = ^ 



g^r, 


(10) 


dS(j) 
dr 


(11) 



at a nonzero, but small, radius r. At the surface, we 
require the p erturbations to be both finite and upwardly 
evanescent. lUnno et al.l ()1989[ ) show that the correct 
boundary conditions are then 

1)^ = 
r 

at the upper boundary of the model. We solve Equations 
(O, ([ni, and ([7]) with the boundary conditions in Equa- 
tions ([HI, dH), (dni), and (HD) using the shooting method 
to obtain all mode periods. 

The middle three panels in Figure [2] display the trans- 
verse displacement eigenfunctions and energy density for 
the n = 1, 5, and 10 modes for £ — 1. The energy den- 
sity illustrates where each mode "lives," that is, what 
portions of the star most affect the mode period. These 
modes live primarily in the core, below the H/He tran- 
sition region, as the energy density declines rapidly in 
the lower pre ssure H layer. This co ntrasts the normal 
ZZ Cetis (Fontaine fc Brassardll2008( ). This is predom- 
inantly because the Brunt-Vaisala frequency in the He 
core is l arger than in the en velope (see Figure [U also 
noted in lAlthaus et al.|[2004[ ). The electron degeneracy 
in all WD cores leads to most of the entropy in the ions, 
yielding N^ ^ A^'^ {kbT / E F){g / Xp) , where A is the ion 
mass and Ep is the electron Fermi energy. Hence, there 
are two reasons why N is relatively larger in an He core 
than in a C/0 core. First, low-mass implies smaller 
Ep oc Af*/"^, and second, 1/A is larger for He than a C/0 
mixture. This shows the power these modes will have in 
probing the core composition once pulsations have been 
detected and accurate periods measured. 

The observability of these modes requires that they be 
driven. One type of driving mechanism requires that 
a portion of the escaping heat flux be converted into 
the mechanical energy of the pulsat fon modes. An ex- 



ample is giv en by t he K-mecha n ism (IDziembowskil 



Dziembowski & Ko esteii 119811 : iDolez fc Vauclaiii 



1977 



1981 



Winget et al.iil982i) . in which the rapid outward increase 



in opacity associated with an ionization zone bottlenecks 
the heat flux. For ZZ Ceti-like pulsations, IBrickhilH 
(il983 i) proposed that the response of the convec tion zone 
itself to the pulsation drives the instability (jBrickhilll 
fl98a[T99lHWu fc Goldreichl[T999l) . For now we consider 
only the convective driving mechanisnfl. If the convec- 
tion zone can thermally adjust on a timescale shorter 
than the pulsation period, V , then the pulsation will be 
damped. This motivates the convective thermal time 
instability criterion V < 87rrth, bcvz, where Tth, bcvz is 
the thermal ti me from the base of the convection zon e 
to the surf ace (BrickhilJ [l99l IWu fc GoldreichI [1991 . 
IWu fc GoTd reich (1999) show (in their Figure 7) this cri- 
terion to be quite accurate in the high T^s (the "blue" 
edge) limit compared to fully nonadiabatic calculations 
on a logg — 8.0 WD. However, the He core WDs are 
at significantly lower gravity, requiring an extrapolation 
in logg. Calculations show that the rapid increase in 
Tth, bcvz due to the deepening convective zone occurs at 

^ The existence of the nuclear burning region inside the mode 
propagation region may additionally drive modes due to the sen- 
sitivity o f the nuclear reactions to temperature, the e-mechanism 
(KawraleilllQSSl V We have yet to investigate this possibility. 
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Fig. 3. — Evolution of stably bur ning low-mass He core s. Evolution proceeds top to bottom. We start our 0.20 Mq model where the 
deviations from the time-dependent ISerenelli et al.l I I2002I) models have been reduced to less than Alogg = 0.2 (see Figure [l}. For our 
0.17, 0.18, and 0.19 Mq models we begin at the maximum Teff turn as we have no models to directly compare. We show our entire 0.15 
and 0.16 Mq models, starting with Mcnv = 1O~'^M0, as our comparison with iPanei et a P |I2007I) show them to be reliable (Figure [l]l. 
For comparison wc exhibit the location of the C/O WDs (the shaded region) and their associated empirical instability strip, the sloped 
dashe d lines (Gia nninas et al. 2007). Candida te objects are drawn from several sources. The bullets are PSR J1911-5958A (Bassa et al] 
I20M 1. PSR J1012+5307 (van Kcrkw iik et al.| [T996: Callanan ct al. 1998), SDSS J0849-I-0445 and SDSS J0822-f2753 (Kilic et al. 2010), 
NLTT 11748 llKawka fc Vermes 2009), an d LP 400-22 (Kawka ct al. 2006). The crosses are those SDSS selected low-mass WDs confirmed 
by MMT spectra from iKilic et al] II2007I ). The circles on the evolution curves represent where the convective thermal time instability 
criterion is met for the 1=1 and n = 1 modes. 



lower Toff for lower logy, nearly 1500 K for logg = 8 to 
logy = 7. The circle points on the logg-Tes evolution 
plots in Figure [3] show where the convective thermal time 
instability criterion is met for i = I and n = 1 modes. 
As is seen, Figure [3] highlights many excellent targets for 
an observational study of He core WDs, as we comment 
in the conclusions. 



4. CONCLUSIONS 

Our work highlights the logy-Toff parame- 
ter space where observable pulsations may be 
present. This reveals at least three pu l sation 

candidates, NLTT 11748 (jKawka fc Venned 120091: 
Steinfadt et all IMol ). SDSS J0822-h2753 (IK ilic et al' 



20101). and PSR J10 12-h5307 (|van Kerkwiik et al. 199 



i 



Callanan et al.|[T998[ ) , all of which should be observed for 
variability on timescales of 200-1000 s. Mode detections 
and a measurement of the mode period spacing would 
provide key evidence for an He core composition and 
large radius, as we predict a period spacing of ~ 90s for 
£ = 1, whereas in normal ZZ Cetis, this num b er is ~ 50 s 
(iKleinman et al.l 119981 : iKanaan et all l2005t iPech et al.l 
|2006[ ). Our results regarding the region of instability 
should be confirmed through future nonadiabatic sta- 
bility analyses, an issue that is beyond the scope of the 
present paper. 

NLTT 11748 is highlighted in Figure E] as a candi- 
date for observ a ble pulsations. Recent observations by 
iSteinfadt et al.l ()201 0l did not find pulsations down to 
a 5 mmag level, however, they did discover it to be 
the first eclipsing He WD system. Given its measured 
logg = 6.20 and Toff = 8540 K, our models predict a to- 
tal mass of 0.17 Mq and envelo pe mass ofS. 15x10"'^ Mp , , 
comparable to that reported bv lKawka fc Venned (|2009l ). 
Our numerical pulsation analysis reveals that for this ob- 
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ject the lowest order 5-mode {i = 1 and n = 1) has a 
period of 245 s but more importantly the mean period 
spacing is 89 s for the £ = 1 modes and 51s for the £ = 2 
modes. Figure [2] clearly shows these modes to preferen- 
tially reside in the core, offering a unique opportunity to 
probe the core composition of an He WD. 

The current candidates were found in surveys 
that target other phenomena: the Sloan Digital 
Sky Surv ey (SDSS; fe ihc ct al. 2007), high veloc- 
ity stars CKa wka et alT l2006t jKawka & Venues' '2009"), 
and companions to pulsars ([van Kcrkwiik ct al. 1996; 
ICahanan et~aLl 119981: fBassa et all I2OO60 . However, the 
survey selection criteria (photometric colors), although 



incomplete, can favor WDs of hig her Teg (Teg > 
11,000K in the SDSS for logg < 6, IKilic et al.l[20oC 
Surveys able to select low-gravity WDs down to Tcs = 
8,000K will significantly impact the study of pulsating 
He core WDs. 
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